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ABSTRACT 

Using the OGLE catalogue of eclipsing binaries, 15 contact binaries were identified 
towards the SMC and the LMC at vertical distances from the Galactic plane between 
300 pc and 10 kpc. Based on the luminosity function calculated for these contact 
binaries, we estimated a frequency of occurrence relative to Main Sequence stars in 
the thick disk at roughly This estimate suffers from the small number statistics, 
but is consistent with the value previously found for the solar neighbourhood. 
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1 INTRODUCTION 



During the past few decades, several attempts have been 
made to determine the local spati al density of con tact bi- 
' naries, with very diversified results. iRucinskil (|2002h (where 
, full references are given and differences discussed) estimated 
■ their local spatial density at (1.02±0.24) x 10~^ pc~^. Later, 
' on the basis of the All Sky Automated Survey (AS AS), a 
[ relative frequency of occurrence (RFO) of one c ontact bi- 
nary among about 500 solar type stars was derived (|Rucinskil 
I2OO6I ). which is in full agreement with the mentioned spa- 
tial density. The RFO for contact binaries at high galactic 
latitudes, however, has been entirely unknown. 



Here we present an estimate for the spatial occurrence 
of contact binary systems relative to Main Sequence (MS) 
stars in the thick disk of our Galaxy, far from the Galactic 
plane. In particular, the luminosity function for contact bi- 
naries in two conical volumes towards the Small Magellanic 
Cloud (SMC) and the Large Magellanic Cloud (LMC) cov- 
ering parts of the thick disk and the halo is determined. The 
results are then compared with the luminosity function for 
MS stars in the same region of the sky. The term "contact 
binaries" is here used as a synonym for W UMa-type eclips- 
ing binaries with orbital periods in a range of 0.22 - 1 days. 
In this paper, for reasons to be explained below, we limit 
ourselves to a subset with periods < 0.45 days. 



2 IDENTIFICATION OF SHORT PERIOD 
CONTACT BINARIES IN THE 
OGLE-CATALOGUE OF ECLIPSING 
BINARIES 

The Optical Gravitational Lensing Experiment (OGLE) was 
intended to detect dark matter in the Milky Way Galaxy us- 
ing the microlensing technique, with the Magellanic Clouds 
and the Galactic Bulge being the main targets of the survey 
jUdalski. Szvmanski and KubiakI Il997l ). As a by-product, 
OGLE provides high quality, long-term photometry that can 
be used to analyse eclipsing bina r y sta rs. Two online cata- 
logues (|Wvrzvkowski et al.ll20o3. I2OO4I ) were used for this 
study. They contain data in the standard photometric BVI 
system for 2850 and 1351 eclipsing binaries, and cover an 
area of 4.6 and 2.4 square degrees of the central parts of the 
LMC and the SMC respectively. The OGLE-II survey has 
a faint hmit of roughly / = 20.5 mag with a corresponding 
error of 0.3 mag; the error becomes smaller for decreasing 
magn itude, reaching a bright li mit of the survey at about 13 
mag (jWvrzvkowski et al.ll2004^ . 

The differentiation between contact binaries and other 
eclipsing binary types in the two OGLE catalogues was car- 
ried out by applying a contact binary criterion based on 
Fo urier analysis o f the light curves, which was introduced 
by IRucinskil (|l997l fl. By means of this criterion and using a 
visual inspection of the remaining light curves, we identified 
10 and 5 contact binaries with orbital periods P < 0.45 d to- 
wards the LMC and the SMC, respectively. The light curves 
are plotted in Figure |31 Furthermore, we used the online 



^ More detailed explanatio ns for this and ot her te chniques used i n 
this paper can be found in l lRucinskilll997h and l lRucinskill200^ . 
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databas e from the MACHO (Massive Compact Halo Object) 
project (|Alcock et al.lll997l ) to confirm the classification of 
the contact binaries in the OGLE catalogue. A direct com- 
parison with the OGLE catalogue was only possible for the 
contact binaries towards the LMC, because there is no MA- 
CHO photometry available for most survey fields covering 
the SMC. The MACHO photometry of the 10 contact bina- 
ries towards the LMC is in good agreement with the results 
obtained from the OGLE survejQ. 

The contact binaries in our sample have orbital periods 
in a range of 0.22 - 0.45 days. T he short period limit is a nat- 
ural cut-off for contact binaries (|Rucinskil2007l ). whereas the 
upper period limit was intentional: On one hand, we wanted 
to be sure of the contact binary classification and avoid semi- 
detached binaries, while on the other hand, we wanted to use 
the simple My = My (log P) calibration and avoid problems 
with uncertain or missing colour indices. Because the final 
RFO estimate is done using absolute magnitude bins of the 
luminosity function, the period limits signify an intentional 
restriction to the low brightness end of the contact binary 
sequence. We will explain the details below. 



1000 




Figure 1. This plot shows the vertical distance from the galac- 
tic mid-plane for the 10 and 5 contact binaries towards the LMC 
(squares) and the SMC (triangles) respectively. The Sun is as- 
sumed to be located 30 pc above the galactic mid-plane. 



3 ABSOLUTE MAGNITUDES & DISTANCES 

Contact binaries show a correlation between the colour (i.e. 
{B — V)o) and th e orbital period P, w hich was first observed 
by EggenI l|l967h . As was shown by iRucinski fc Duerbeckl 
this correlation can be used to determine an ab- 
solute magnitude calibration My = Mv {log P,{B — V)o)', 
an even simple r version, My = My (log P), is applicable 
l|Rucinskill2006l ) when a restriction to short periods is added. 
In this research, the absolute magnitude calibration My = 
-1.5(±0.8) - 12.0(±2.0)logP was used, which is applica- 
ble for P < 0.562 d (|Rucinskii,200a '). While applying this 
calibration to the contact binary sample, we arbitrarily as- 
sumed an uncertainty in absolute magnitude of 0.3 mag; the 
formal uncertainties in the calibration were larger and prob- 
ably overestimated. Using this absolute magnitude calibra- 
tion and the relation between the distance and the apparent 



V-Mv + 5-A^ 



the dis- 



and absolute magnitudes, d = dex( 
tance from the sun and the galactic height for the 15 contact 
binaries were estimated. The extinction towards the LMC 
and S MC was obtained from (jSchleeel. Finkbeiner &: Davis! 
1 19981 ): we assumed Ay — 3.1Eb-v- Figure [1] illustrates 
the vertical distance from the galactic mid-plane to the 15 
contact binaries as a function of the absolute magnitude, 
whereas Table [1] lists the numerical results. 



4 THE AMPLITUDE DISTRIBUTION AND 
THE LUMINOSITY FUNCTION FOR THE 
15 CONTACT BINARIES 

Small-amplitude systems often remain undetected in vari- 
able star searches when the photometric error is too large 
for detection of a small photometric variability. The number 



^ As an exception, the OGLE photometry in V-band for the con- 
tact binary OGLE050905. 22-693315.1 was found to have a gross 
error of roughly 3 magnitudes. 



of missed small amplitude systems was estimated by com- 
paring our sample with a theore tical model of t he amplitude 
distribution for contact binaries (|Rucinskil200ll ) . This model 
predicts the amplitude distribution as a function of the de- 
gree of contact / and the mass ratio of the two companions 
q = ^ for a sample of contact binaries. In this research, we 
assumed / = 0.25, what appears to be the favoured degree 
of contact, and a fiat d istribution of q. We applied this model 
in the same way as in (i Rucinskillioog ). where more informa- 
tion can be found. Following this approach, we multiplied 
the number of contact binaries by factors of 1.6 and 1.5 for 
the LMC and the SMC respectively; both factors have an 
uncertainty of about 20%. 

In order to determine the luminosity function for the 
15 contact binaries, we simply used 1-mag wide absolute- 
magnitude bins by dividing the (corrected) number of stars 
in each bin by the respective volume V = ^{r2 — ^1)0. C 
denotes the fraction of the sphere covered by the survey and 
is equal to 1.09 x lO""* for the LMC and 5.82 x 10"^ for the 
SMC. The results are given in table [2] 



5 STELLAR DENSITY MODEL AND 

LUMINOSITY FUNCTION FOR MS STARS 

In order to establish the relative numbers of contact binaries, 
in relation to MS stars, we attempted to derive the predicted 
numbers of stars in the same search volumes as those of 
the OGLE and MACHO surveys. A ccord ing to the pioneer- 
ing work of iGilmore fc ReidI l| 19831 ) and iBahcall fc Soneiral 
( 1984 ), the stellar distribution in the Milky Way can be 
modelled by a double exponential thin and thick disk and 
a spheroidal halo. Here, we describe the stellar density by 
such a three component model 

n{z,R) = nthin{z,R) + nthick{z,R) +nsp{z,R), 

where {z,R) are the galacto-centric coordinates. The thick 
and the thin disks are modelled by 



nthin{z, R) 



no • (1 — Cthick ■ 
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Table 1. The short period contact binaries towards the LMC and the SMC. The table lists the name of the star (OGLE 
convention), the Magellanic Cloud field, the orbital period, the apparent magnitude in V-band, the absolute magnitude My, 
the distance along the line of sight and the vertical distance above the Galactic mid-plane. 

name field period (d) V (mag) My {mag) distance (fcpc) vertical distance (fcpc) 



OGLE003835 


24-735413 


2 


(SI) 


SMC 





26909 


15.89 


5.3 


1.22 ± 0.17 





82 


± 





12 


OGLE005431 


85-723510 


9 


(S2) 


SMC 





33125 


16.44 


4.3 


2.59 ± 0.36 


1 


77 


± 





26 


OGLE005846 


47-724315 


3 


{S3) 


SMC 





33244 


15.61 


4.2 


1.78 ± 0.25 


1 


21 


± 





18 


OGLE004619 


65-725056 


2 


(34) 


SMC 





37663 


13.55 


3.6 


0.93 ± 0.13 





62 


± 





09 


OGLE004849 


85-725554 


8 


(35) 


SMC 





43328 


18.67 


2.9 


13.78 ±2.01 


9 


54 


± 


1 


44 


OGLE054003 


85-703837 


3 


(LI) 


LMC 





24091 


15.48 


5.9 


0.73 ± 0.10 





37 


± 





06 


OGLE050542 


01-691725 


9 


(L2) 


LMC 





27755 


16.31 


5.2 


1.51 ± 0.21 





79 


± 





12 


OGLE051932 


28-694633 


4 


{L3) 


LMC 





27871 


16.64 


5.2 


1.78 ± 0.25 





94 


± 





14 


OGLE053540 


37-695413 


9 


{L4) 


LMC 





29159 


16.30 


4.9 


1.69 ± 0.23 





89 


± 





14 


OGLE053251 


73-700256 


2 


(L5) 


LMC 





29718 


17.71 


4.8 


3.39 ± 0.47 


1 


82 


± 





28 


OGLE053916 


98-700903 


3 


(L6) 


LMC 





31132 


15.35 


4.6 


1.28 ± 0.18 





67 


± 





10 


OGLE050905 


22-693315 


1 


{L7) 


LMC 





32563 


17.00 


4.3 


3.05 ± 0.42 


1 


63 


± 





25 


OGLE053247 


54-694403 


1 


{L8) 


LMC 





33104 


17.34 


4.3 


3.71 ± 0.51 


1 


99 


± 





30 


OGLE053539 


86-694759 


5 


(L9) 


LMC 





36844 


15.09 


3.7 


1.70 ±0.24 





90 


± 





14 


OGLE054701. 


27-705623.3 (LIO) 


LMC 





42700 


18.18 


2.9 


10.01 ± 1.41 


5 


44 


± 





82 



Table 2. The luminosity function for the 15 close contact bina- 
ries towards the LMC and the SMC. The limiting distances are 
determined from d = dex((V — My + 5 — Ay)/b) with apparent 
magnitude limits V G (13±0.3, 21±0.3). The errors in the number 
of stars n are calculated according to the Poisson statistics, which 
contributes to the uncertainty in cji as the main factor. Further- 
more, we arbitrarily set the number of stars equal to one for the 
SMC bin centred at My = 6 with zero star counts; an assumption 
which is consistent with the 1 — a Poisson error. 



Cloud 



My 





{mag) 


{kpc) 


{kpc) 




/ stars \ 


SMC 


3 


1.19 


30.00 


1.5 


0.23 ±0.21 




4 


0.75 


18.93 


4.5 


2.71 ± 1.71 




5 


0.48 


11.94 


1.5 


3.62 ±3.31 




6 


0.30 


7.54 


1.5 


14.39 ± 13.18 


LMC 


3 
4 
5 
6 


1.13 
0.71 
0.45 
0.28 


28.41 
17.93 
11.31 
7.14 


1.6 
4.8 
8.0 
1.6 


0.15 ±0.14 
1.82 ± 1.12 
12.10 ±6.59 
9.63 ±8.60 



X exp 



Zthi 



nthick{z,R) = no • Cthic 



X exp 



Zthick 



exp 



exp 



hthin 

R — -Ro 



hthick 



where Zthin and Zthick denote the scale heights and hthin 
and hthick the scale lengths for the thin and the thick disk 
respectively. Furthermore, no is the local stellar density, and 
Cthick and Cap are the local-density normalizations of the 
thick disk and the halo relative to the thin disk. In this 
research, we use a projected de Vaucouleurs spheroid (jYoungj 
ll976l : lBahcall fc Soneiralll984l ) of the form 



nsp{R') = no ■ Csp 
X exp 



10.093 ( ^ ) + 10.093 



Here, R' is given by 



where k is the axis ratio of the de Vaucouleurs spheroid. 

In our model, we used Zthin ~ 300 ± 50 pc, zthick ~ 
1000 ± 200 pc, htMn = 2.5 ± 1 kpc, hthick = 3.5 ± 
1 kpc, Cthick = 0.05 ± 0.03, Csp = 0.0015 ± 0.001, k = 
0.55 ± 0.1 and _Ro = 8 ± 0.5 kpc. These jaarameters were 
chosen in consideration of the work of Gilmore & ReidI 
1983), Kuiikcn & Gilmore (1989), Fux & Martinet (199jx 



[5 



Bieiiavme fc SechaudI (Il997l ). iLarsen fc Humphreys! (120031 ) 
and lJuric et al.l i 2005h . The rather large uncertainties in the 
parameters reflect ranges in the results obtained in the men- 
tioned p apers. 

Wie len. Jahreiss fc Kriiged (|l983l ) derived a widely ac- 
cepted luminosity function for MS stars in the solar neigh- 
bourhood using a sample of nearby stars, which we used 
to determine the stellar density in the solar neighbourhood 
no as a function of absolute magnitude. With those results 
and the stellar distribution n{z, R) as described above, we 
calculated the mean stellar density of MS stars in each vol- 
ume corresponding to the four, 1-mag wide, absolute magni- 
tude bins. By doing so, we were able to compare the derived 
contact binary luminosity function with the expected MS 
luminosity function for the same portion of the sky. The 
luminosity functions towards the LMC and the SMC are 
plotted in figures [2] and [3] respectively. The contact binary 
luminosity functions were found to be best approximated by 
the MS luminosity function, if the latter was divided by 650 
and 600 for the LMC and the SMC, respectively. 



6 CONCLUSIONS 

We conclude that each of the two samples yields the relative 
frequency of one contact binary among about 600 MS stars 
in the two conical volumes towards the SMC and the LMC 
at galactic latitudes of —44° and —33° respectively. These 
estimates have uncertainties of roughly 50% due to the large 
uncertainties arising from small number statistics and to the 
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Figure 2. Luminosity function towards the LMC. Tlie continuous 
line represents the contact binary luminosity function, whereas 
the dashed line represents the MS luminosity function divided by 
a factor of 650. 



^ 1.00 
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My-bin (mag) 



Figure 3. Luminosity function towards the SMC. The same as 
in figure [2l except for the scaling factor of 600 instead of 650. 



still uncertain parameters of the spheroidal component of 
the Galaxy. This component has been explicitly accounted 
for in the previous section and is estimated to contribute to 
the total number of stars in the SMC and the LMC search 
volumes at levels of about 30% and 20%, respectively. The 
contribution varies with the Mv- Specifically, for the SMC, it 
changes from 38% for the M^-bin with the largest distances 
to 17% for the M^-bin centred at 6 mag; for the LMC, it 
changes from 29% to 10% in the same range of Mv- 

The results on the relative frequency of contact bi- 
nary stars at large galacto-centric distances is consistent 
with estimates for the solar neighbourhood in previous work 
l|Rucinskil 120061 ) ■ These are the very first data available on 
the contact binary distribution at high galactic latitudes. 
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Figure 4. The light curves of the 15 close contact binaries. The fitting was done according to the approach of l|Rucinskilll993^ (i.e. using 
only the first 5 cosine and the first sine Fourier-coefficients). 



